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We discuss different operator definitions of the transverse-momentum dependent (TMD) parton 
densities from the point of view of their renormalization-group (RG) properties and UV evolution. 
We also consider the structure of the gauge links (Wilson lines) in these operator definitions and 
examine the role of the soft factor in the factorization formula within the TMD approach to semi- 
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inclusive processes. 

PACS numbers: ll.10.Jj, 12.38.Bx, 13.60.Hb, 13.87.Fh 



The candidates for the golden measurement at the EIC are the spin-dependent Sivers function fyj,, as well as the 
' unpolarized quark distribution f\. The proposed silver candidates are the transversity, the Boer-Mulders, and the 
Collins functions. All these objects are transverse-momentum dependent (TMD) parton densities (termed in the 
following for simplicity, TMDs) that describe the inner structure of hadrons by taking into account the partonic 
pH : longitudinal and also the transverse degrees of freedom [H-Q. 

i— r The QCD factorization formula for a semi-inclusive structure function is expected to have the following symbolic 
form (modulo power corrections) 

>• F~H®F D ®F F ®S. (1) 
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This expression contains the hard part H, the distribution and fragmentation transverse-momentum dependent func- 
tions Td and J~f, respectively, and the soft part S, which is introduced in order to take care of light-cone singularities 
• , that cannot be cured by regularization. While the hard part H can be evaluated order by order within the framework 
of perturbative QCD, the TMDs are essentially nonperturbative entities and have to be modeled or extracted from 
the data. Thus, little can be said about these objects, except their general operator definition and their evolution 
behavior with respect to the energy and rapidity scales. Contrary to the fully collinear, i.e., integrated PDFs, the 
• • ' explicit operator definition of TMDs is quite problematic, owing to the appearance of extra (light-cone) singularities, 
. £h ' anc ^ remams under active investigation (see, e.g., (§-[r7|). 

Let us concentrate on the following two definitions 1 of an unpolarized quark TMD 

h(x, k ± |n) = \ J ^^e- ik+r+i ^ (hMC,Z±)[C,Z±;°c-,z ± }t 

x[oo~,£_L;oo~,oo_ L ]j' 7+ [oo~, oo ± ; oo~, Oji[oo~, ± ; 0~, J„V>i(Cr, 0±)\h) (2) 



and its counter-part 



fi(x, kx|u)=/i(a:, kjn-^i;) . (3) 



following the so-called Amsterdam notations, e.g., [19j . while other TMDs can be analyzed within the same framework. 
Gauge invariance is ensured by means of path-ordered contour-dependent Wilson-line operators (gauge links) with the 
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1 A generalized definition, which includes into the Wilson lines the spin-dependent Pauli term -F M1/ [7 M , y u ], is worked out in Ref. [l8j | — see 
also 03]. 



2 



generic form [y, x\T] = "Pexp — ig JJ™ dz^iz)]- where A = t a A a , and one has to distinguish between longitudinal 

[ j ][n.v] and transverse [ , ]ri] gauge links 0, On the other hand, the definition used in lattice simulations 

contains, in contrast to the (semi-)infmite Wilson lines in Eqs. @ and ([3]), the finite, i.e., the direct gauge link joining 
the two quark fields [U^ll]. 

In the tree-approximation, the "distribution of a quark in a quark" is normalized as 

/f oc (x, k ± \n) = S(l-x)S^(k ± ) 

and the integration over kj_ yields formally the usual collinear (integrated) PDF 

Jd 2 k ± A (Z, k x |n) - f$ Iee (x) =6(l-x). 

However, definition ([2]l. taken literally, produces — beyond the tree-level — certain pathological divergences, which 
belong to one of the following three classes: 

1. Usual UV- singularities ~ i from integrations over loop momenta, which can be removed by using the standard 
R— operation. 

2. Pure rapidity divergences, which only appear in the unintegrated case. They cancel in the integrated distributions, 
but they are present in the TMD case giving rise to logarithmic and double-logarithmic terms of the form 
~ In C, In 2 C; they have to be resummed by a consistent procedure. 

3. Overlapping divergences, which contain both UV and soft singularities simultaneously: ~ ^ In meaning that 
the UV pole e _1 mixes with a "soft" divergence, regularized by the auxiliary parameter £. 

In addition, one may encounter combinations of the above types of singularities ensuing from the soft factor. Apart 
from the UV divergences, all other singularities originate from uncompensated light-cone artifacts, which stem either 
from the lightlike gauge links (in covariant gauges), or from specific terms in the gluon propagator in (singular) light- 
cone axial gauges. While the singularities of the second class may be simply regularized by a rapidity cutoff — which 
is "separated" from other variables — the effect of the third class is more severe, because these singularities affect the 
UV-renormalization procedure, change the anomalous dimensions, and modify, therefore, the RG-evolution. 

In order to avoid the above-mentioned problems, the following approaches have been proposed in the literature: 

1. Shift in covariant gauges the gauge links off the light-cone: v 2 ^ 0, or use instead the non- lightlike axial gauge 
(v ■ A) — 5]. This amounts to definition ([3]), but it may cause problems in properly deriving factorization 11 



2. Stay on the light-cone, Eq. @, but subtract some specific soft factor R, which is defined in such a way as 
to exactly cancel the extra divergences (23l - [25j : Eq. @ is substituted by the "subtracted" function fi(n) — > 
h{n)-R-\ 

3. Perform a direct regularization of the light-cone singularities in the gluon propagator [26!] l/q + — > ^/[q + ](r]) , 
where rj is an additional dimensional parameter [12j. In this case, a generalized renormalization is in order, 
which is formally equivalent to multiplying the TMD PDF by a particular soft factor [27| . 



4. Use the light-cone axial gauge, but supply it with the Mandelstam-Leibbrandt (ML) pole prescription [28l l29j: 
l/q + — > l/{q + + iQq~) or equivalently q~/(q + q~ + i0). Now the overlapping singularities do not appear 
from the outset — at least at the level of the one-loop order — while the contribution of the soft factor is reduced 
to unity, rendering the gauge- invariant definition valid [3] . 

Let us present the UV evolution equations for the above definitions. The off-the-light-cone TMD does not 
contain overlapping singularities. Therefore, the renormalization-group equation reads 0, l30| 

Mx" fi(x,k±,n\v) = 7lc fi(x,k±, (i\v) , 7lc = 7 + 0(a 2 s ) , (4) 

d\i A TT 

where 7lc is the anomalous dimension of the bilocal quark operator in the light-cone gauge. If one factorizes out the 
soft contribution R v , as proposed in Ref. Q, then the anomalous dimension changes and one has 

h(v,fi) -> fi(v,[i) ■ R^ 1 , fi-^ [f%(v,fi) ■ Ry 1 ] = (7lc ~ 7r) fi{ v Tf) ' Ry 1 i ( 5 ) 
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where 7r is the one- loop anomalous dimension of the soft factor R v . 

In contrast, the anomalous dimension of the "light-cone" TMD, before subtraction, deviates from 7lc and this 
deviation is determined by the cusp anomalous dimension (T^ |: 



V— h( n ) = (7lc ~ 7cus P ) fi(n) . (6) 



A 

d[i 

Hence, the generalized rcnormalization procedure restores the "broken" anomalous dimensions, so that one finds 



/i(n) -> h(n) ■ R- 1 , fiA- [h(n) ■ R- 1 ] = 7LC [h{n) ■ R' 1 ] . (7) 

In the light-cone gauge with the Mandelstam-Leibbrandt prescription, Eq. ^ supplemented by the soft factor R 
yields ab initio an anomalous dimension without lightlike artifacts: 

^ [/l( " )ML ' R ^ = ^ A(n)ML = 7LC [/l(n)ML ' R " 1] = 7LC A(n)ML ' (8) 

Thus, evolution equations can be established, while the evolution with respect to the rapidity variable — either £, or r\ 
(depending on the approach applied) — constitutes a separate task. 

To conclude, let us sketch a couple of important but still un(re)solved problems. 

(i) Factorization of semi-inclusive processes: an all-order factorization (in a covariant gauge) was studied in Ref. 
0, but definition © was used which contains off-the-light-cone longitudinal gauge links. An explicit proof of a 
factorization theorem with lightlike longitudinal gauge links in the TMD PDFs is to our knowledge still lacking. 

(ii) Relationship between TMDs and collinear PDFs: the generalized definition with lightlikc longitudinal gauge 
links ([7]) does indeed yield after integration an a;— dependent distribution function that obeys the DGLAP evolution 
equation 

d 2 k±_ fi(x,k±,n I n) = fi(x) , M-y- fi(x,fJ-) = /Cdglap ® fi(%, A*) ■ (9) 

dfi 

The reason is that the overlapping singularities in the real and the virtual gluon contributions cancel against each 
other after carrying out the k^-integration. In contrast, for an axial but not-lightlike gauge, one gets a function 
containing longitudinal gauge links off-the-light-cone along the vector v = (v + , v~ , Ox)j i-c, 

d 2 kj_ fi(x,k±,/j,\v) = fi(x,n\v) , (i — fi(x,(i\v) = K v ® f\{x,n\v) , K v ^ /Cdglap ■ (10) 

The RG-properties of this object differ from those of the TMD PDF with lightlike longitudinal gauge links (which 
fulfills the DGLAP equation shown above). This difference leads, in particular, to different evolution equations and 
different phenomenological implications. 
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